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We compare standard perturbation theory with the polaron transformation for 
non-linear transport of electrons through a two-level system. For weak electron- 
phonon coupling and large bias, there is good agreement between both approaches. 
This regime has recently been explored in experiments in double quantum dots. 



1 Introduction 

Semiconductor coupled quantum dots are well-defined artificial few-particle sys- 
tems for the study of electron-electron correlations, quantum coherence, and quan- 
tum dissipation. These effects become visible in non-linear electron transport on 
an p,eV energy scale when parameters like the interdot coupling and the coupling 
to external leads can be controlled by e.g., external gate-voltages.lil 

Double quantum dots in the regime of strong Coulomb blockade can be tuned 
into a regime that is governed by a (pseudo) spin-boson modelcl (dissipative two- 
level system) 



H= £ -a z ~ ' 
2 

Q Q 



T c o- x + ^a z A + ^2uj Q a f Q aQ, A :=^ ffQ (a~Q + Gq) , (1) 



where one additional 'transport' electron tunnels between a left (L) and a right (R) 
dot with energy difference e and inter-dot coupling T c , where a z — \L) (L\ — \R) (R\ 
and a x — \L)(R\ + |i£)(£|. Here, ujq are the frequencies of phonons, and the gQ 
denote interaction constants. Although not exactly solvable, the model is quite 
well understood for closed systemsS (isolated dots with one additional electron). 
The coupling to external leads offers the possibility to study its non-equilibrium 
properties, such as the inelastic stationary current through the dots. 

2 Equations of Motion 

We describe the dynamics of the double dot by a reduced statistical operator p(t), 
allowing for an additional 'empty' state and tunneling from a left reservoir at rate 
Tt, into the left dot, and from the right dot to the,right reservoir at rate T^. Lowest 
order perturbation theory in these rates yieldsEra 

d 

-Q t PLL{t) = -lT c [pLR\tj - PRL\t)\ + T £ [1 - PLh{t) - PRR\t)\ 

d 

g^PRR(t) = ~%T C [pRL \t) - pLR{t)] - r R p R R(t). 
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For the remaining equation for the off-diagonal element plr. = P*rli one nas ^° 
choose between perturbation theory in gQ (weak coupling, PER), or in T c in a 
polaron-transformed frame (strong coupling, POL). In general, no exact solution 
of the model is available: this is the case even for coupling to one bosonic mode 
only {gQ oc 5q,q , Rabi Hamiltonian). 

For the spin-boson problem with Tr/l — 0, it is well-known that POL is 
equivalent to a double-path integral 'non-interacting blip approximation' (NIBA) 
that works well for zero bias e = but for e ^ does not coincide with PER at small 
couplings and low temperatures. In the following, we compare both approaches for 
Fr/l 7^ and find nearly perfect agreement for very large e » T c , a regime that 
has been tested experimentally recently!!! 

The standard Born and Markov approximation with respect to A yields 

J t PlTit) = [*'e - 7„ - Ik/2] PLR (t) + \%T C - «J_] p RR (t) - [zT c - 6+] p LL (t). 
Here, the rates are 

> := 2n^p(A) coth (/3A/2) , S ± := -ff ~p(A) coth ((3A/2) T §|p(A) 
p(Lu):=J2\9Q\ 2 S("-"QY (2) 



Q 



where A := y/e 2 + 4T C 2 is the energy difference of the hybridized levels, and f3 = 
1/ksT the inverse phonon equilibrium bath temperature. Note that beside the 
off-diagonal decoherence rate 7 P , there appear terms oc 5± in the diagonals which 
below turn out to be important for the stationary current. 

On the other hand, the polaron transformationa leads to an integral cquati 



tion 



V -fC{t~t') PLR {t') 



P P L° R L {t) = - dt'e^S) 
Jo 

+ iT c {C(t - t') PLL (t') - C*(t - t') PRR (t')} 



where 



3 Stationary Current 



C(t) := cxp <{ - / duj^- [(1 - coswt)coth(/3w/2) +isinujt} \ . (3) 



If an electron tunnels between two lateral dots, the interaction with 3d piezoelectric 
acoustic phonons leads to an orthogonality catastrophe ('bcaon shake up effect') 
which is determined by an effective phonon density of statescl 



p(uj) = guj I 1 - -2 sin j j e-»^9(cj), (4) 

showing oscillations on a scale LOd '•= c s /d, where c s is the speed of sound, d the 
distance between the centers of the two dots, g the dimensionless coupling constant, 



mbllproc: submitted to World Scientific on February 1, 2008 2 



and lo c a high-energy cut-off. We argue that the assumption of sharply localized 
positions between which the additional electron tunnels is justified by the strong 
intra-dot electron -electron repulsion. 

We compare results for the stationary electron current 7 stat = — e2T c lmp*L R (z — 
0) as obtained from Laplace transforming the equations of motion within both PER 
and POL for small elect ron-phonon coupling. The result 

j = -eT L T R G+ 

stat T L G- + {T L +T R )G+-T L T R () 

has to be used with either 

q(per) := 2TJm _JIk^_ (6) 
is - 7 P - T R /2 

or 

G< POL) := 2T c Im , G (POL) := 2T e Im , (7) 

+ i + (i/2)r R c e ' - l + (l/2)r R C £ ' y) 

where C £ := J °° dte l£t C(t). 

The function C' e is obtained for arbitrary large coupling constant g by a double 
numerical integration. For our purpose here, it is sufficient to expand C £ as 

G £ = -+ 7 (e) + 0(. g 2 ), 7(^:=~^[coth(/?| E |/2)+sgn( £ )]. (8) 
4 Discussion 

In both PER and POL the current Eq. (|) is an infinite sum of contributions 
G+ (= istat/e in lowest order in T c ) and G_ and therefore to infinite order in the 
inter-dot coupling T c . Note that both expressions for G± coincide for vanishing 
electron-phonon coupling g. 

We first point out that PER works in the correct eigenstate base of the hy- 
bridized system (level splitting A), whereas the energy scale e in POL is that of 
the two isolated dots (T c = 0). This difference reflects the general dilemma of 
two-level-boson Hamiltonians: either one is in the correct base of the hybridized 
two-level system and perturbative in g, or one starts from the 'shifted oscillator' po- 
laron picture that becomes correct for T c = 0. In fact, the polaron (NIBA) approach 
does not coincide with standard damping theoryu because it does not incorporate 
the square-root, non-perturbative in T c hybridization form of A = e 2 + 4T C 2 . 

However, for large \e\ 3> T c , A — > |e|, and POL and PER should coincide. This 
is indeed the case for small g where 

2 T R /2 + 1 {±e) 
-G ± -2T C ry4 + £2 . (9) 

Thus, for large |e|, the polaron approach works well even down to very low tem- 
peratures and small coupling constants. In particular, in the spontaneous emission 
regime of large positive e the agreement is very good. This regime is most in- 
teresting at. low temperatures T and has been tested in the non-linear transport 
experimental in great detail. 
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In Fig. 1, wc compare both approaches and show the calculated stationary 
current, using realistic parameters. In the vicinity of e = 0, the agreement between 
POL and PER gets worse as could be expected from the above argument. At higher 
temperatures, the agreement gets slightly better on the absorption side e < 0. We 
conclude that our findings for the 'open' spin-boson model are in agreement with 
standard spin-boson physics. Furthermore, we see that for large bias |e| ^> T c both 
perturbation theory and the polaron transformation approach practically coincide 
for small coupling g. 

This work was supported by the DFG, Project BR 1528/4, and the EPSRC, 
GR/R44690/01. 
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